Spin conservation and Fermi liquid near a ferromagnetic quantum critical point 
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We propose a new low-energy theory for itinerant fermions near a ferromagnetic quantum critical 
point. We show that the full low-energy model includes, in addition to conventional interaction 
via spin fluctuations, another type of interaction, whose presence is crucial for the theory to satisfy 
5*17(2) spin conservation. We demonstrate the consistency between a loop- wise expansion and 
a Fermi liquid description for the full model. We further show that, prior to the ferromagnetic 
instability, the system develops a Pomeranchuk-type instability into a state with zero magnetization 
but with p-wave deformations of the Fermi surfaces of spin-up and -down electrons (a spin nematic). 

PACS numbers: 71.10. Ay, 71.10 Pm 



Introduction Conservation laws are powerful tools to 
study the behavior of interacting fermions. The Hamilto- 
nian for a system of fermions with non-rclativistic inter- 
actions is SU (2) invariant and preserves the total charge 
and spin of particles. In this Letter, we analyze the conse- 
quences of spin conservation for effective low-energy the- 
ories, which describe the behavior of itinerant systems 
near a q = Pomeranchuk instability. In particular, we 
focus on the effective theory of a ferromagnetic quantum 
critical point (FM QCP). It has been long assumed that 
a FM QCP is well described by the spin-fermion model 
— (SFM) with interaction only in the spin channel, medi- 
ated by small q collective spin excitations (for a review, 
see Refs.ii). 

Here, we report two results. First, we show that the 
SFM is fundamentally incomplete as it violates spin con- 
servation. We obtain a proper low-energy model which 
obeys spin conservation and show that in this new model 
the interaction is mediated by both spin and charge 
bosons; the charge component of the interaction is of the 
same strength as the spin component. We show that the 
perturbation theory for this model is consistent with the 
Fermi liquid (FL) description, which is based on conser- 
vation laws. This consistency does not hold in the SFM, 
as it has been recently emphasized in Ref. 0. 

Second, we show that, within the new model, the sys- 
tem develops an instability toward a p— wave spin ne- 
matic state,— 1^ in which the net magnetization is absent 
but the Fermi surfaces of spin-up and spin-down electrons 
are shifted in opposite directions. For not too long-range 
interaction in the spin channel, this instability occurs not 
only before the Stoner instability into a ferromagnetic 
state but also before other known instabilities near a FM 
QCP -a first order transition, helical magnetic order, or 
p-wave superconductivity. ^ii^i^iiS 

A p— wave spin nematic state, first introduced by 
Hirsch,'^ has been studied in detail in recent years.- It 
was identified as the ground state of a particular lattice 
model,— but whether such a state can be realized in a ro- 
tationally and SU (2)-spin invariant model was not clear. 



Our result shows that such a state emerges in an ordinary 
isotropic FL tuned to a FM QCP. 

Spin conservation in the SFM. The SFM model was 
introduced as a minimal model near a magnetic QCPf^iii 
It describes low-energy fermions with an effective spin- 
spin interaction 
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Here, q = k — p, = cj^ — Wp are the relative momenta 
and Matsubara frequencies, respectively, S measures the 
distance from the QCP, v is the density of states, and Z 
is the fermionic residue. The superscript Q implies that 
Eq. ([1]) is obtained by integrating out fermions outside 
the Fermi surface (FS) but contains no contribution from 
fermions at the FS - the latter is obtained within the FL 
theory which takes P^ an input. The main structure 
of Eq. HJ) is obtained already within the RPA^ d^d"^ (see 
Fig. [T]i); an additional factor oi 1/Z arises from renor- 
malization of beyond the RPApi^ 

The effects of the interaction P^ can be analyzed in two 
ways - either within a loop- wise expansion in P^ or with 
the FL theory built upon conservation laws. In a model 
which satisfies conservation laws, the results of the two 
approaches must be identical. We first show that they 
are not in the SFM. For definiteness, we consider below 
the case of D = 2, but our results are also valid for D = 
3. Loop- wise expansions of bosonic and fermionic self- 
energies have been performed by a number of authors, 
and the results relevant to our discussion arei^^ 
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Observe that the static uniform spin susceptibility is not 
renormalized by the interaction at low energies. 

To obtain the results for m* /m and Xs in the FL the- 
ory, one has to convert P^^ into the Landau function 
5a,3;^5(k, p) = Z'^{m* /m)T^p.^^{'k, 0; p, 0), obtain charge 
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FIG. 1: a) Ladder diagrams for F . Summing up these 
diagrams for a contact interaction U (dashed line) and re- 
arranging spin indices, one obtains the effective spin-spin in- 
teraction near the critical point (wavy line), where Um/2TT « 
1. The rest of the diagrams (not shown) eliminate the insta- 
bility in the charge channel, b) Renormalization of the spin- 
fermion vertex. The last two diagrams are the Aslamazov- 
Larkin-type (AL) contributions which restore spin conserva- 
tion, c) The full effective interaction p'^'^"" with contributions 
from the direct spin-spin interaction and AL-type processes. 




the particles on the FS (|k| = |p| — kp, uj^ — — 0), 
and use famihar FL relations m*/m = 1 + gc,i and 
Xs = X°(l + ffc,i)/(l + gs,o)- Carrying out this proce- 
dure, we find gc,i — A and gs.o = —A/3, i.e.. 



The result for m*/m is consistent with Eq. ((2]) but the 
result for Xs is not. Moreover, Eq. ^ shows that Xs 
has a spurious divergence at A = 3, i.e., at finite distance 
from the QCP. 

This disagreement between the two approaches is re- 
lated to a violation of spin conservation in the SFM. To 
see this, we compute the spin and charge susceptibili- 
ties in the SFM at a small but finite frequency and 
zero momentum q. The Ward identities require that 



Xc,s{q = 0,n) = 0. We show that Xs{q = 0,n) 0, 
although Xc{q = 0, ri) = 0. 

The effect of the interaction with collective modes on 
the spin and charge susceptibilities can be parameterized 
as 

Xc,.(g = 0,l])=xL(l-^), (4) 

where Ac,s are the fully rcnormalizcd fermion-boson ver- 
tices normalized to Acs = 1 for non-interacting fermions. 
Self-consistency requires that Ac,s are obtained within 
the same RPA approximation as Eq. ^ itself. The 
RPA series [cf. Fig. [2))] is geometric; the n-th order 
term is given by (A/(l -I- A))" for the charge vertex and 
(-1/3)"(A/(1 + A))" for the spin vertexJI The difference 
comes from spin summation: the charge vertex contains 
a combination 5ai3<^a-t ■ f?(5/3 = 3(5^^, while the spin vertex 
contains crl^pcfa-y-ffsp = — f^^- Summing up the series, we 
find that Ac = 1 + A and A^ = (l-HA)/(l-h4A/3). Substi- 
tuting Acs into Eq. (|4]), we find that charge susceptibility 
Xsil — 0,a;) vanishes, as it should, but Xsiq — 0, does 
not, which implies that spin conservation is violated. 

This violation has a physical explanatiouf^i In the 
SFM, electron spins are split into spins of itinerant elec- 
trons, s = c^CTq/3C^ and spins S of collective bosons. In 
the SU (2) symmetric case, the fermion-boson coupling 
s • S flips s^. As a result, s is not conserved separately 
from S. 

Because of spin non-conservation, the SFM was ar- 
guedii to be applicable only in the regime when conser- 
vation laws do not matter, i.e., when relevant momenta 
are much larger than relevant frequencies. Some physical 
properties described by the SFM, e.g., mass renormaliza- 
ton, superconducting and other instabilities of a param- 
agnetic state, come from fermions from that range and 
are therefore adequately described within the SFMi^ii^ 
Still, the violation of spin conservation implies that the 
SFM model is fundamentally incomplete as the SU{2) 
invariance is not broken in a paramagnetic phase at any 
energy. A correct low-energy model, which includes all 
relevant interactions, should obey spin conservation. 

Spin- conserving low- energy model. We now con- 
struct such a model by going beyond the ladder approxi- 
mation for the vertex corrections Ac and As. It turns out 
that the important diagrams are the Aslamazov-Larkin 
(AL) type-diagrams encountered, e.g., in the theory of 
superconducting fluctuations^^ To second order, there 
are two such diagrams shown in the last line of Fig. [1] 
6). For the charge vertex, the two diagrams cancel each 
other For the spin vertex, they add up and contribute 
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Here, G(k, = Zl\iijj — efe(m/m*)] is the Green's func- 
tion of FL quasiparticles, = vp(^ — kp), m* /m = 
= 1 + A, 6lim = Z(l,m), and f(q,17) = 

(5 + {aqY + \^\/vFq) ■ Although /al is formally a 
two- loop correction in F, its contribution is actually of 
the same order as from the one-loop correction. To see 
this, we observe that relevant ep in Eq. ([5]) are small (of 
order Og) while q is finite. Hence, the integration over 
ep involves efi^ectively only the first two Green's functions 
in the first line. For the poles of these two Green's func- 
tions to be in different half-planes, the range of uj must 
be bounded by the external frequency — fio < < 
Hence, one can neglect uj and fio in the remainder of 
the integral. The angular integrations now factorize and 
each of them gives a factor of [m* /m)/vpq. The re- 
maining integrals over q and have a scaling property 
d{n/vFq)f'^{q,^/vFq) ^ J dqT{q,0) - A. A more 
accurate calculation yields /al = (4/3)A/(l + A) (the 
denominator comes from power-counting factors of Z). 
Combining /al with the one-loop ladder vertex correc- 
tion, which contributes — (1/3)A/(1 + A), we find that 
the effective one-loop correction to becomes A/(l-|-A), 
which is exactly the same result as the correction to Ac. It 
can be readily shown that the effective nth other contri- 
bution to As now is (A/ (1+A))", again the same as for Ac. 



Summing up the series, we obtain that now As = 1 + A, 
which means that Xs(0, fio — ^ 0) = and thus spin con- 
servation is restored. Other two-loop vertex corrections 
are small as max{r2FL//?F, ^^o/^^fl}, where ilpL oc S^^'^ 
is the upper boundary for the FL description. 

For completeness, we also analyzed spin conservation 
right at a FM QCP, where the fermionic self-energy 
S(w) oc w^/^. The calculations at QCP are more in- 
volved, and to verify the spin Ward identity one has to 
solve an integral equation for the full vertex. For the spin 
vertex Tg we find the same result as above: the AL-type 
diagrams restore spin conservation. 

Once we have established that AL blocks play the role 
of an effective single bosonic propagator, it is obvious 
that these blocks also give additional contributions to 
the FL vertex r2^.^^(k, 0; p, 0) and, hence, to Landau 
parameters. These contributions, shown Fig.[T]c), add an 
extra term of the structure {3dasSp-y^(TaS'<^i3'y)F{\i+p) to 
r^^.^5(k, 0; p, 0). Although F(k -I- p) is a rather complex 
function of its argument, one can show that its partial 
harmonics are the same as those of the RPA-like kernel 
l/{6 + a^(k -I- p)^). This statement is true for harmonics 
with n ^ A which are the only ones we need. Combining 
Eq. Il]) with the AL contribution, we obtain for the full 
vertex near a FM QCP 
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Two comments are now in order. First, the AL term 
in Eq. ([6]) has both spin and charge components and the 
latter is by no means negligible. In other words, a spin- 
conserving low-energy model contains singular interac- 
tions in both spin and charge channels. Second, the AL 
term has the same structure as the direct spin-fluctuation 
exchange term, but becomes singular at small total mo- 
mentum (k w — p) rather than at small transferred mo- 
mentum (k w p). 



Fermi liquid near a FM QCP. 



We 



construct the new Landau function ff^iL^Ck, p) 



Z\m*/m)T'^^^^}^,ik, 0; p, 0) = g, 
from the effective interaction F^'^"'' 
dard spin algebra, we obtain 
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where 6 — Z(k, p). Evaluating g^^'f and 5s"o\ we obtain 
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A for A ^ 1, hence 



m*/m = A; Xs = X. 
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These two results are now in full agreement with those 
of the loop- wise expansion, Eq. 

The FL theory with the Landau function from Eq. ([7]) 
is rather non-trivial. First, not only 5c,i but all 
charge components with n <^ \ diverge near critical- 
ity in the same way: g^^^l — (3/27r) d9 cos{n9) / {S -\- 
4a'^kl,sm'^ 9/2) « A (see Fig. The divergence of 

an infinite set of Landau parameters also exists near a 
charge QCP,i4 This does not lead to dramatic conse- 
quences as the divergence of gc,n>i is compensated by 
that of m* /m = 1 -|- gc i, and the charge susceptibilities 
Xc,n = Xluim*/m)/{l+gif) = x° „(1+A)/(1+A) = 
remain intact. 

Second, and more important, we see from Eq. ([7]) 
that all odd spin Landau parameters are negative and 
diverge at the FM QCP. Since a Pomeranchuk insta- 
bility occurs when the corresponding Landau parameter 
approaches —1, it implies that a spin Pomeranchuk in- 
stability with finite angular momentum occurs before the 
isotropic, Stoner-like instability. Indeed, to leading order 
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FIG. 2: (Color on-line.) Schematic evolution of the partial 
components of the full Landau function [Eq. ^] vs the in- 
verse mass renormalization parameter A~^. All charge and 
even spin Landau components are positive but diverge near 
a FM QCP, while odd spin components are negative. The 
spin nematic component gs,i reaches the Pomeranchuk criti- 
cal value of —1 first. 



full 
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in 1/A, we have from Eq. ^ g^"^^ = A(5„^2m - f A^„,2m+i- 
For n = 2m + 1, g*"^' — —5 A/3, and gl^^ approaches —1 
at A = 3/5. To discriminate between states with differ- 
ent angular momenta, one needs to evaluate the Landau 
parameters beyond the leading order in 1/A. Doing so, 
we find, quite naturally, that g^"}' approaches —1 first. 
Therefore, the leading instability corresponds to a p- wave 
spin nematic state. Such a state in itinerant Fermi sys- 
tems has recently attracted significant attention.^ Two 
different phases (a and /3) -the analogs of A and B phases 
of 3He- have been identified, but only the a phase is actu- 
ally stable. This phase is characterized by a p— wave spin 
order parameter = I]k /('^)("^kT'^kT - Ck^Ckj,) cos6'k, 
where f{k) is sharply peaked at the FS. This state has 
zero magnetization but the centers of FSs of spin- up and 
spin-down electrons are separated by finite momentum. 

As we have already mentioned, a p— wave instabil- 
ity competes with more "conventional" instabilities of 
a FL near the FM QCP - p— wave superconductivity. 



a first order transition and a transition into a helical 
magnetic phase. Which one occurs first depends on 
the parameters of the model. Controllable calculations 
are possible only when the range of interaction in the 
spin channel is large, i.e., for akp ^ 1. According 
to Ref. i, the most relevant out of "conventional" in- 
stabilities is the first order transition, which occurs at 
6 — (5^,. ~ 0.21/{akp). The p-wave Pomeranchuk insta- 
bility occurs at (5 = S^,. = 25/16{akF)^- Although S^^ is 
smaller than Sl^ by power-counting, the two critical val- 
ues coincide at large enough value of akp « 7.44, which 
implies that for a realistic case oi akp ^ 1, the p-wave 
instability comes first. 

Finally, we notice that SU (2) spin invariance is only 
approximate in real systems and broken by, e.g., spin- 
orbit (SO) coupling. This breaks the exact relation be- 
tween the ladder and AL-type vertex corrections, but the 
AL term in F^ still remains of the same order as the bare 
spin interaction, and, as long as SO coupling is small, the 
system still undergoes a p— wave spin-nematic instability 
at A - 1. 

To conclude, we considered the behavior of an SU{2) 
invariant itinerant fermionic system near a ferromagnetic 
quantum critical point. We showed that the spin-fermion 
model, traditionally used to describe this behavior, does 
not satisfy SU (2) spin conservation. We obtained a spin- 
and charge-conserving low-energy model for a FM QCP. 
This model has singular interactions in both spin and 
charge channels, as opposed to the SFM which considers 
interaction only in the spin channel. We constructed a 
FL theory for such a model and found that the system 
undergoes a Pomeranchuk instability into a p— wave spin 
nematic state before a FM QCP is reached. 
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